A consistent quasiparticle picture of Quark-Gluon plasma and the velocity of sound 
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velocity of sound in a hot quark-gluon plasma. The results are compared with the recent lattice 
data. 
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The study of the Equation of State (EOS) of strongly interacting matter is of great importance. It can play a vital 
role for a deeper understanding of many unresolved areas of QCD physics, such as the study of quark-gluon plasma 
(QGP), quark confinement, phase structure of QCD at high temperature and finite chemical potential, to name a few. 

The study of QGP has some formidable problems. Due to the convergence problem in the low energy domain, the 
widely used perturbative QCD calculation fails badly to give any insight. The most accepted method in this area is the 
lattice QCD calculation, though the method works better with light dynamical quarks. Moreover, there are difficulties 
with finite baryon chemical potentials. On the other hand, various phenomenological models of thermodynamical and 
hydrodynamical types 0,llllllllEllllll[ll[lJG2G3 are used as they are easier to handle compared to the 
lattice. In all these models, a quasiparticle description is used with a background mean field. Here it is worth noting 
that a quasiparticle description is expected to be applicable, as long as the spectral functions for quarks and gluons 
resemble qualitatively the asymptotic forms found from HTL (Hard Thermal Loop) perturbative calculations [14( . 

The calculations in the quasiparticle description usually start with a Hamiltonian. It was pointed out by Biro et al. 
|15| that in a Hamiltonian approach, due to the dependence of the Hamiltonian on the thermodynamical characteristics 
of the surrounding matter, there may be some thermodynamical inconsistency within the method which has to be 
rectified. To get rid of these inconsistencies, one has to incorporate certain conditions in the quasiparticle model. 

In all the studies of the EOS, the sound velocity (it), which relates the pressure and the energy as P = u 2 E, plays an 
important role. Specifically, the QGP is expected to approach the canonical value of u 2 — 1/3, at high temperatures 
and/or densities. An understanding of the quantitative behaviour of the approach to asymptotic freedom is very 
important, as it is related to the problem of quark confinement. 

Another important quantity related to the study of QGP is the quark number susceptibility (QNS), which, being 
related to charge fluctuation, is of direct experimental relevance. The importance of QNS lies in the fact that it can 
be used as an independent check on the theoretical models which try to explain the lattice resul ts Qa. 

There are numerous results of quark number susceptibility available from lattice calculations (l6lll7lll8l |. Also, There 
have been recent reports of lattice studies where the velocity of sound in a pure gluon plasma at finite temperature 
has been calculated fl9| . Our motivation here is to understand what kind of quasiparticle picture can reproduce these 
results. 

In the present work, we have studied the behaviour of quark number susceptibility and the speed of sound as the 
quark-gluon plasma moves from the non-perturbative to the perturbative region of QCD. The density dependent 
quark mass (DDQM) model [20l l2lj. along with the condition for thermodynamical consistency, has been used for 
the present study. 

The DDQM model aims to explain confinement from a phenomenological point of view. The MIT bag model [l| is 
so far the simplest option available to explain the confinement phenomenon, though it has some serious shortcomings. 
The DDQM model is a simple but effective alternate approach to confinement. The main idea behind the DDQM 
model is to make the mass of the "free" quarks infinitely heavy (confinement) while at small distances, effective quark 
mass should be small or zero j2^|. This behaviour is similar to what is known as "Archimedes Effect" [23t l24j . 

In the DDQM model, confinement is achieved through a parametrization of the effective quark mass as a function of 
density. In the low density limit, all the interactions are taken care of by the mass, so that it becomes infinitely large 
to prevent the occurrence of a single quark in the asymptotic state 20] . The thermodynamical quantities related to 
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the quark sector at finite temperature are calculated from the postulated mass relation and a self-consistent number 
density relation. At finite temperature the quark mass is defined as |2l| . 

m{T,n q ) = + , (f) 

where nf(n~) is the number density of the i-th flavour quark (antiquark) and n q is the net quark number density. 
These number densities arc to be calculated self-consistcntly from the given relation: 

nf = — / dpp 2 {eMT-\^T^)] + l}-\ (2) 
n Jo 

with, 

e l (p,n q ,T) = [p 2 +m 2 (T,n q ) + A. l (p)} 1 / 2 (3) 

where Aj(p) takes care of the divergent term present in the mass at finite temperature. In the zero density limit, 
Aj(0) = as the mass takes care of all interactions. The expression for A$(p) is given by [2l| : 

dAiiPi) = \— x WHB^mw* + dPi (4) 

with Cj = rii/n q . The energy density and the pressure of the quark system are then readily calculated from 

o poo 

S «=J2— dp^e^n^aexpiT-^e.-^ + iy' + lexpiT-^e.+^ + l}- 1 )} (5) 
i— ^ •'0 



P i = o E — / Po-^PfiexpfT-^ei - Mi)] + + {exp^fe + Mi)] + I}" 1 ) 



(6) 



In the gluonic sector, energy density and pressure have the form 



45 47T 

Here, a c is the effective gluon-gluon coupling. The form of a c is given by plj : 



54B7T 3 

(< q^T„t) 2 



" c = — 02 ^ „\2 ( 8 ) 



where, 



+ ni)<Q 2 >i +n g < Q 2 > g 
<Q 2 >n g ,T=y ^— (9) 



< Q 2 >i and < Q 2 > g are the thermal average of the three momentum squared of quarks of flavour i and gluons. 

< Q 2 > g corresponds to non-interacting gas of gluons, n g being the gluon number density [2lll25|. 
So the total energy density and pressure of the quark-gluon plasma is, 

Sqgp = £ q + S g , Pqgp =P q + P B , (10) 

In the regions of low temperature and density, the energy density and pressure of the gluonic sector becomes negative 
which can be interpreted as the signal of the formation of gluonic condensates. As long as the gluons remain in the 
condensed phase they do not contribute explicitly to the thermodynamic quantities (they are not active degrees of 
freedom of the system). So, £ g = V g =0 for a c > jg plj . 
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The DDQM model has been widely used in the literature for the study of dense quark matter |26|, 1271 l28l I29L 
|33,|^|, the thermodynamical inconsistency notwithstanding Let us now modify the above EOS incorporating 
the conditions of thermodynamical consistency. The prescription for ensuring thermodynamical consistency is to 
introduce a "background" field $ which should satisfy the following conditions 



gtn. 

13; 



j=u,d 

where (i, is the degeneracy factor and v k j is the average quasiparticle occupation number of j-th flavour quark with 
momentum k: 

v kj = ({explT-^ej - ^)] + l}" 1 + {exp^- 1 ^ + ^)] + l}" 1 )] (13) 
In the present case, the above two conditions become, 

ar+E4V = ' (14) 

^ + E^ = 0- (16) 

(s) 

where is the scalar quasiparticle density of j-th flavour quark, 

nf = ^^-2 (16) 

In the present paper we have considered a two flavour QGP system. Since u and and d quarks have equal masses, our 
system, in essence, consists of one type of quasiparticles only (except for the degeneracy factor 2 for two flavours). 
So, the quasiparticle quark mass depends on temperature and quasiparticle density through the scalar density only 
and the expression for the background field <& reduces to 

* = -/ E n j )dn i ( 17 ) 

j=u,d 

Carrying out the integration and identifying the integration constant with $ = at T=0, jib = 0, background field $ 
can be obtained. At zero temperature, <E>, naturally, depends on the baryon density (jib) alone. With this background 
field, we define thermodynamically consistent energy density and pressure as, 

Sqgp =S q + £ g + $(n B ,T) (18) 
P Q GP = P q + P g -Hn B iT) (19) 

At zero temperature and density, the total energy density of the QGP becomes B, the energy of an empty bag. On 
the other hand, the pressure of QGP becomes zero and we do not encounter the unphysical negative pressure as found 
in most of the MIT bag like phenomenological models. It is to be noted here that there is no bag pressure in this 
model as such. Since the mass of the quark becomes infinitely large in the density going to zero limit, a negative bag 
pressure would amount to double counting |20| . 

We define the baryon (rij = o) and isospin (n, = 3) densities and the corresponding susceptibilities as: 

_ T dlnZ _ T d 2 lnZ 

where, Z is the partition function and thermodynamic potential is defined as il = —Pqgp = —TlnZ. Also 
Mi=o = M« + and [i i=i = 

In figure 1 and 2 we have plotted y.t/yf ft as a function of T/T c . Here X3 is the iso- vector susceptibility and Xfft 
is the free field value; for details see [32, |33| . It is to be noted here that the value of the critical temperature in 
figure 1 and 2 are different. In figure 1 the T c is taken from the consideration of the commonly held belief that the 
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QCD critical temperature is between 150-200 MeV; in our calculation it is taken to be 200 MeV and represented by 
dash-dot curve. In figure 2 the critical temperature has been calculated from the DDQM (taking the zero pressure 
limit, beyond which quark matter system will become unstable) and shown as continuous line in the graph. The points 
are from lattice calculation 0] for different values of the valence quark mass. As seen from the plot, our result is in 
qualitative agreement with the lattice result in the region near the critical temperature. But at higher temperature it 
approaches the ideal gas limit faster than what we get from the lattice calculation. Here it is worth mentioning that 
in our calculation current quark mass is zero and the dynamic quark mass varies with temperature. 

The variation of $ with temperature at three different baryon densities for B 1 / 4 =145 MeV is shown in figure 3. 
The value of $ is maximum for ub = and it decreases for increasing density. This behaviour of of <& could have been 
anticipated, as the system is expected to be closer to the perturbative domain at higher densities. Near the transition 
temperature T c , $ shows a rapid variation and then saturates at high temperature. In fact, around the saturation 
region the effect of 4> on thermodynamic quantities becomes negligibly small. In figure 4, we have plotted (£-3p)/T 4 
with temperature. The two curves are for the cases with and without the background field B 1 / 4 and ub being fixed 
at 145 MeV and 3no (where no = 0.17 fm~ 3 is the nuclear matter saturation density) respectively. The continuous and 
the broken lines correspond to the cases with and without the background field contribution, respectively. Inclusion 
of the background field into the model increases the critical temperature by 29% in the DDQM. The non-perturbative 
effects are dominant near the transition temperature. The inclusion of thermodynamic consistency certainly enhances 
these effects, as seen in figure 4. Similar behaviour is obtained for the other values of B as well. 

Our main motivation in this paper has been to study the velocity of sound, which is given by (|f )s- But in the 
present paper we have plotted (§f )tj as it is easier to compute in the present model. Moreover, it was pointed out 
in rcf.[2l| that the difference between (§f )t an d (§f )s does not exceed 10%. The figures 5, 6, 7 and 8 show the 
variation (§f )t as a function of temperature for different densities at B 1/,4 =145 MeV. We have also presented the 
recent results from lattice calculation [l^. As mentioned earlier, the dash-dot curves in figures 6 and 8 are plotted 
with T c equal to 200 MeV. The continuous and the broken lines in figures 5 and 7 correspond to the cases with and 
without the background field contribution, respectively, the T c being calculated from the DDQM. 

As seen from all the plots, (§f )t is greatly modified by the contribution from the background field around T c ; $ 
has a dominant effect near T c . Moreover, as <E> goes over to B at zero temperature and density, the change in B affects 
the sound velocity much more than that due to the change in density. It is to be noted here that the lattice results for 
sound velocity [li| presented here are for gluonic contribution only, whereas our results are for quark gluon plasma. 
The discontinuity in the plot is due to the built-in second order phase transition (gluon condensate) in the present 
EOS. 

Here we would like to point out that in figures 7 and 8, we have compared our results for the speed of sound in 
QGP at non-zero ub with the lattice values at = 0. These comparisons indicate that the behaviour near T c is 
governed mainly by the physics of confinement. In our model, the presence of quarks (along with the gluons) has very 
small effect on the in-medium behaviour of sound velocity. 

To conclude, we have studied the velocity of sound in the QGP with a thermodynamically consistent DDQM 
model. This is a one parameter model and the dependence on the parameter is also seen to be very weak. This is 
the remarkable feature of this model. Incorporation of thermodynamic consistency produces qualitatively and semi- 
quantitatively, the features of the lattice data and also explains all the necessary features of the non-perturbative 
to perturbative transition of the QCD physics. Our study shows that the incorporation of consistency condition is 
essential to understand the nonperturbative behaviour near transition temperature. The result obtained for QNS falls 
within the range of lattice data. The result for the sound velocity of quark gluon plasma for zero baryonic density at 
critical temperature 200 MeV is in qualitative agreement with the lattice result for pure gluon plasma. Our result of 
sound- velocity of QGP for the finite baryonic density may be considered as a prediction to be confirmed by the future 
lattice calculations. 
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FIG. 1: Variation of Xs/xfft with temperature for B 1 ^ 4 =145 MeV ,ns=0 and T c = 200MeV. Dash-dot curve represents 
our model calculation. Different symbols like circle, down-triangle, up-triangle, diamond and the square represent lattice data 
corresponding to the valence quark mass 1.0, 0.75, 0.5, 0.3, 0.1, in units of T c , respectively Hfj. 




FIG. 2: Variation of Xs/Xfft with temperature for B 1//4 =145 MeV ,ns=0 and T c =330 MeV, calculated from the model. 
Solid curve represents our model calculation. Different symbols like circle, down-triangle, up-triangle, diamond and the square 
represent lattice data corresponding to the valence quark mass 1.0, 0.75, 0.5, 0.3, 0.1, in units of T c , respectively 0. 




FIG. 3: Variation of $ with temperature for B 1 ^ 4 =145 MeV. Dash-dot, solid and dotted curves are for densities ns=0, 
ns=3no and ns=4no respectively, no being the nuclear matter saturation density. Critical temperature for riB—O is 330 MeV, 
for n s =3n it is 290 MeV and for n B = 4n it is 260 MeV. 
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FIG. 6: Variation of (fj)r with temperature for B 1//4 =145 MeV, ns=0 and T c = 200MeV '. Lattice data (circles) are taken 
from (l9|. Dash-dot line is for our model. 



0.3 - 



3 0.25 



0.2 - 



0.15 



1 1 1 1 1 1 1 1 


i 1 i 1 






- 


— 


/ 

/ / 

k / 

A/ / 
* 

1 
i 


- 


- JO 








6 




o 1 


— 


1 i , i , i 





4 

T/T 



6 



7 



FIG. 8: Variation of (§j)t with temperature for B 1/ ^ 4 = 145 MeV, n_B=3no and T c = 200MeV. Lattice data (circles) are taken 
from [l9|. Dash-dot and the broken lines correspond to the cases with and without <fi contributions, respectively. 



